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ABSTRACT 

This work focus on the first results on the identification and characteriza- 
tion of periodic plasma density fiuctuations in the outer corona, observed in 
STEREO-A CORl white-light image time series. A 2D reconstruction of the 
spatial distribution and temporal evolution of the coronal fiuctuation power has 
been performed over the whole plane of the sky, from 1.4 to 4.0 Rq. The adopted 
diagnostic tool is based on wavelet transforms. This technique, with respect to 
the standard Fourier analysis, has the advantage of localizing non-persistent fiuc- 
tuating features and exploring the variations of the relating wavelet power in both 
space and time. The map of the variance of the coronal brightness clearly outlines 
intermittent, spatially coherent fiuctuating features, localized along, or adjacent 
to, the strongest magnetic field lines. In most cases they do not correspond to 
the coronal structures visible in the brightness maps. The results obtained pro- 
vide a scenario in which the solar corona shows quasi-periodic, non-stationary 
density variations characterized by a wide range of temporal and spatial scales 
and strongly confined by the magnetic field topology. In addition, structures 
fiuctuating with larger power are larger in size and evolve more slowly. The char- 
acteristic periodicities of the fiuctuations are comparable to their lifetimes. This 
suggests either that plasma fiuctuations lasting only one or two wave periods 
and initially characterized by a single dominant periodicity, rapidly decay into a 
turbulent mixed fiow via nonlinear interactions with other plasma modes, or that 
they are damped by thermal conduction. The periodic non-stationary coronal 
fluctuations outlined by the closed field lines at low and mid latitudes might be 
associated with the existence of slow standing magneto-acoustic waves excited 
by the convective supergranular motion. The fiuctuating ray-like structures ob- 
served along open field lines appear to be linked either to the intermittent nature 
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of the processes underlying the generation of magnetic reconnection in the polar 
regions or to oscillatory transverse displacements of the coronal ray itself. 

Subject headings: methods: data analysis — magnetohydrodynamics (MHD) — 
Sun: corona — Sun: magnetic topology — Sun: oscillations — waves 



Introduction 



MagnetoHydroDynamic (MHD) oscillations and waves in the coronal plasma are of great 
interest in solar physics due to the cru 
heating and solar wind acceleration (e.g 



Cranmer et al. 


1999; 


Walsh & Ireland 


2003; 


Ofman 



20041). Coronal waves are also irnportant as a probe for MHD coronal seismology (e.g. 
Uchidal Il970t [Roberts et al.l Il984j : iNakariakov fc Verwichtd |2005| ) : the comparison of the 
observed properties of these phenomena, such as periods, wavelengths, amplitudes, spatial 
and temporal scales and time evolution, with theoretical models may lead to cast light 
on physical quantities of the solar corona, such as the magnetic field and the transport 
coefficients, not yet accurately known from the observational point of view. 

Periodic and quasi-periodic phenomena in the solar corona w ere first observed with 
ground-based radio instrumentation i n the 80s (I AschwandenI 119871 ). The launch of SOHO 
(SOlar and Heliospheric Observator y. iDomingo et al.lll995[) . in 1995, and TRACE (Transi- 
tion Region And Coronal Explorer, Handy et al.l 19991 ). in 1998, opened a new era in the 
investigation of coronal waves and transients. The highly improved spatial and temporal 
resolution of the new instrumentation carried by these spacecraft allowed the identification, 
almost ubiquitously in the corona, of a large variety of MHD waves: standing fast kink distur- 
bances, fast sausage and slow longitudinal acoustic modes, propagating slow and fast waves, 
and so on. They permeate individual structures of the inner corona, such as loops, polar 
plumes and helmet streamers, with different spatial and temporal scales and periods rang- 
ing from a few seconds to several minutes. Longer period (8 — 27 h) oscillations have been 



observed with EIT/SOHO (Extreme ultraviolet Imagi ng Telescope, Delaboudiniere et al. 



19951 ) in the inner corona in a large-scale EUV filament (jFoullon et al.ll2004j ). UVCS/SOHO 
(Ultraviolet Coronagraph Spectrometer, iKohl et al.lll995l ) has allowed to discover density 
fluctuations in the outer corona, in the region where the solar wind is accelerated. Slow 
magneto-acoustic waves with periods of a fe w minutes (6 - - IQm in) have been detected i n 
a coronal hole between 1.5 Rq and 2.1 Rq by lOfman et al.l (120001 ) and [Morgan et al.l ( 120041 ). 
Large-scale density fluctuations with longer periods (from a few ho urs to a few days) have 



been discovered in both the fast wind at high latitude at 2.1 Rq ( iBemporad et al. 



Telloni et al.ll2009bl ) and in the slow wind at mid and low latitudes at 1.1 Rq (iTelloni et al. 



2008 
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2009al ). The density fluctuations found in tlie slow wind can be interpreted as due to tlie 
existence of inliomogeneities locally generated by nonlinear interactions and carried out by 
the expanding coronal plasma. In alternative they can be also ascribed to outward propagat- 
ing magneto-acoustic waves characterized by wavelengths compatible with the supergranular 
scale in the convection zone. 

In spite of the wealth of the existing observations, the investigation on coronal plasma 
fluctuations is still fragmentary and affected by heavy limitations. Up to now, only the 
dominant temporal scales of periodic and quasi-periodic variability have been detected, and 
in general information on their temporal evolution is still lacking. Moreover, due to the 
limited instantaneous Field of View (FOV) of the spectrometers used in the previous studies, 
only a small fraction of the Plane Of the Sky (POS) has been investigated at a given time. 
This did not allow us to derive the spatial distribution of the plasma fluctuations. 

Aim of this paper is to overcome these major limitations in the study of the periodic 
and quasi-periodic variability of th e outer corona, by exploiting the improved observational 



capab ilities of STEREO-A CORl flThompson et al.ll2003l : Howard et al.ll2008l : iKaiser et al. 



20081 ). launched in 2006, and by adopting a diagnostic technique apt to provide both spa- 
tial and temporal information on coronal plasma fluctuations. This is achieved by analyzing 
long-term time series of the total brightness derived from the high-cadence white-light images 
obtained with CORl and by adopting the wavelet transforms as analysis technique, which 
allow us to separate spatial and temporal variations. The detailed study of the spatiotem- 
poral evolution of the non-stationary plasma fluctuations, observed at different timescales, 
allows a 2D reconstruction of the wavelet power which maps the outer corona from 1.4 to 
4.0 Rq in a systematic manner and over a time interval of 9 days. This approach allows us 
to investigate the correlation between coronal periodic and quasi-periodic phenomena and 
magnetic structures and to attempt an identification of the origin, excitation and dissipation 
mechanisms of the detected variability. 

The layout of the paper is as follows: description of the observations (§ 2), description 
of the analysis of the data and presentation of the results (§ 3), discussion (§ 4) and inter- 
pretation (§ 5) of the results, conclusions (§ 6). The appendix reports methodological details 
of the wavelet technique. 



2. Observations 



The present analysis is based on white-light images of the outer solar corona obtained 
during the minimum of solar activity in 2008 with the coronagraph CORl-A, which is 
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part of the SECCHI (Sun-Earth Connection Coronal and Hehospheric Investigation) instru- 
ment suite aboard the 'ahead' spacecraft of the twin STEREO (Solar Terrestrial RElations 
Observatory) satellites. CORl-A observes the solar corona from 1.4 to 4.0 Rq, providing 
1024x1024 pixel images with a spatial resolution of 7.5 arcsec (~ 0.008 i?©) per pixel. 

The entire set of observations consists of = 2592 exposures acquired with a cadence 
of 5 mm, in the 9 days interval between April 14, 2008 at 00:05 UT and April 22, 2008 at 
23:55 UT. During the selected observational run, the Sun was in an extremely quiet period 
with very few sunspots, small active regions and no evidence of eruptive phenomena. 

Continuous temporal series of white-light total brightness are derived for each pixel of 
the CORl-A images by using the IDL routine SECCHLPREP in the SolarSoft library. This 
routine allows also the calibration of the data and the subtraction of the background. The 
total brightness is computed in mean solar brightness Bq units by combining the intensity 
radiation / of a rapid polarization se quence triplet measured at polarizer positions 0, 120, 



240 de grees: B — |(/o + -^120 + -^240) (lBillingslll966l ). As an example. Fig. [T| shows the total 



brightness of the solar corona on Ap ril 18, 2008 at 21:50 UT; the SOHO/MDI (Michelson 
Doppler Imager. IScherrer et al.lll995l ) magnetogram (as seen from the STEREO-A satellite 
on April 19, 2008 at 00:04 UT) is superimposed to the white-light image. In the same figure, 
the magnetic configuration of the corona is ou tlined by extrapolating the coronal magnetic 
field lines according to the model developed by ISchrijver fc DeRosal (120031 ) (employed in the 
SolarSoft package named PFSS); this model is based on the evolving full-Sun Carrington 
maps of the photospheric magnetic field measured by the SOHO/MDI instrument. The 
coronal magnetic field is extrapolated from the photospheric one via the Potential-Field 
Source-Surface (PFSS) approximation, in which the field is assumed potential in the coronal 
volume between the photosphere and a spherical source surface located at a height of 2.5 Rq. 
Since the coronal field models are provided at six-hour cadence by the online database of 
PFSS, Fig. [1] shows the magnetic configuration closest in time to the total brightness frame 
of April 18, 2008 at 21:50 UT. The extrapolations have been generated taking into account 
the instantaneous line of sight of STEREO-A by considering the Carrington latitude and 
longitude values at each time. 



3. Analysis and results 

In this paper we adopt the wavelet analysis with the aim to identify periodic and quasi- 
periodic density fluctuations over the whole outer solar corona and to explore their fine spatial 
structure and temporal evolution at different timescales. In recent years, this technique, 
applied to observations of the inner corona, led to the detection of periodic phenomena 
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Fig. 1. — Total brightness of the outer solar corona, from 1.4 to 4.0 Rq, derived from the 
CORl-A observations obtained on April 18, 2008 at 21:50 UT; the extrapolated magnetic 
field lines and the SOHO/MDI disk magnetogram, obtained on April 19, 2008 at 00:04 UT 
as seen from STEREO-A are also shown; the white dashed line indicates the edge of the 
STEREO occulter. 
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such as: i) outward propagating; o scillations in an active region bright loop with periods of 



3 — 7 min (iDe Moortel et a. 



quiet c orona ([Morgan et al 



region (IBloomfield et al.l 120041 ) . 



20001) . ii) sporadic oscillations with periods of 7 — 8 min in the 



20041 ) and iii) coherent fluctuations in a weak magnetic network 



The wavelet analysis is in fact a powerful tool for analyzing non-stationary time series or 
time series with expected localized variations of power. By decomposing a time series into the 
time-frequency space, it is possible to determine both the do minant f req uencies /p erio d icit ies 
of the variability and the related power modulation in time ( Torrence fc Compdll998l ). The 
main benefit of the wavelet over the Fourier analysis resides in the capability of identifying 
the periodic components of non-stationary signals. By combining the results obtained at 
multiple locations, the spatial distribution and temporal variability of a given data set can be 
reconstructed. Details on the wavelet formalism and methodology, with particular attention 
to the application to the present analysis, are reported in the Appendix. 

We applied the wavelet transform to the 9-days time series of the total brightness derived 
for each pixel of the CORl-A images. Trends affecting the temporal series in this interval 
have been previously removed; the effect of the solar rotation is filtered out by applying a 
Fourier highpass band filter with a c utoff frequency of 2.8 x 10~^ H z, by taking into account 



the coronal rotation rate derived by iGiordano fc Mancusd (|2008[ ) . The long uninterrupted 
time series and the short temporal sampling, 5 mm, allows the detection of periodicities 
ranging from 10 min (the Nyquist-Shannon period) up to 10 h, the cutoff timescale of the 
highpass filter. 

Since the Sun was very quiet during the observational interval, without any evidence of 
eruptive or transient phenomena, the data are smoothly varying without abrupt jumps or 
sharp peaks. As an example, the left panel of Fig. [2^ shows the total brightness measured at 
2.05 i?0 and 212° counterclockwise from the North pole (which refers to a ray-like structure, 
see § 3.1) before (black curve) and after (red curve) the filtering process. The right panel of 
Fig. [2^ shows the result of a Gaussian smoothing of the red curve which acts as a lowpass 
filter for period smaller than 2.5 h and allows us to cast a preliminary rough glance at the 
total brightness variability behavior. Highly unstable fluctuations whose amplitudes vary 
even significantly with time (strong impulses alternate with phases of damped amplitude) 
are present all along the observational run. This behavior is noticed in arch-shaped structures 
as well (see § 3.1). Figure [2)d shows the wavelet power spectrum derived from the filtered 
brightness time series reported as a red curve in the left panel of Fig. [2^, by using the Paul 
mother wavelet (see the appendix). The ordinate is the Fourier period r, corresponding to 
the scale s of the wavelet, the abscissa the time of observation. The contour plot indicates the 
power of the total brightness fluctuations derived by means of the wavelet transforms, that 
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is the variance of the total brightness time series detected by STEREO- A CORl. The cross- 
hatched area below the Cone Of Influence (COI) line (continuous line in Fig. ^) is the region 
of the wavelet spectrum where edge effects, due to the finite-length of the time series, are 
significant. The shorter periodicities found above the COI line are due to the non-stationary 
periodic signals present in the data. The thick dashed contours delimit regions where the 
wavelet power is significant, that is they identify real periodic sign als at a confidence level of 
99%, as determined via a 2 degree of freedom distribution (see iTorrence fc Compdll998l . 
and references therein as well as the appendix). 

As shown in Fig. [2]d, most of the power is concentrated in the 4 — 8 /i periodicity range. 
In this band its value exhibits significant changes, with a peak on April 20, 2008 at about 
12:00 UT. As expected, negligible power is found at periods larger than ten hours, where 
frequency contributions have been filtered out. Hence, we are confident that the high-power 
features identified during the observational run are not to be ascribed, for instance, to the 
passage of coronal features across the instantaneous FOV during the solar rotation. 

A vertical slice through the wavelet plot of Fig. is a measure of the spectrum at a 
given instant. The average of the local wavelet power spectra over the observational interval 
gives the global wavelet spectrum (as defined in Eq. IA4p shown in Fig. [2t. It is worth 
noting that the global wavelet spectrum approaches th e results that can be obtained with 



the Fourier spectrum of the time series (jPercivall 119951 ). The global spectrum is above the 



99% confidence level (red dotted line in Fig. in correspondence with the periodicities in 
the 4 — 8h band where the highest power is found. 

In order to examine how the wavelet power change throughout the observational run, 
the power spectrum has been divided in four different periodicity bands (0.5 — 1 h, 1 — 2h, 
2 — Ah and 4 — 8 /i). For each band the average wavelet power, that is the average of the 
power over all the periodicities within a selected range [si - S2], has been computed (Eq. 
IA5|) . This quantity is a measure of the average variance of the total brightness as a function 
of time. The total brightness power exhibits significant modulations with timescales that are 
of the same order of the periodicities considered: the primary peaks in each band are one 
or two periods wide. Moreover, the longer the periods considered the slower the temporal 
evolution of the total brightness variance (Fig. |2J1). We will discuss a possible physical 
interpretation of this temporally intermittent behavior in § 4 and § 5. 

Since the results of the wavelet analysis is strongly dependent on the choice of the 
wavelet function (see the appendix), the present results (Fig. [2)d) are compared with those 
derived using the Morlet wavelet (Fig. Wp)- 



The same features appear to be significant in both wavelet power spectra, approximately 
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Fig. 2.— (a) left: Coronal total brightness detected on April 18, 2008 at 2.05 i?© and 212° 
counterclockwise from the North pole (black curve); the red curve has been obtained by 
filtering out any trend, including that due to the solar rotation, from the original time series 
(black curve), (a) right: Gaussian smoothing of the filtered time series (red curve on the 
left panel) which acts as lowpass filter for periods smaller than 2.5 h; note the finer scale 
resolution of the ordinate axis, (b) Wavelet power spectrum of the filtered time series, using 
the Paul wavelet function; the ordinate is the Fourier period r in h; the abscissa is the time of 
observation; the continuous line is the Cone Of Influence (COI) and delimits a cross-hatched 
area where the power spectrum is not significant, whilst the thick dashed contours delimit 
areas where the wavelet power is significant with a 99% statistical confidence, for a white- 
noise process, (c) Wavelet power spectrum averaged over the entire observational interval 
for the dataset reported in (a) (i.e. global wavelet spectrum); the red dashed line is the 99% 
statistical confidence for a white-noise Drocess. (d) Averages of the wavelet Dower over the 
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Fig. 3. — Same as Fig. |2l in the case of the analysis performed with the Morlet wavelet. 
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at the same time-periodicity locations and with comparable power. However, since the Paul 
mother wavelet is more sharply defined in time than the Morlet function (Fig. [7]), the peaks 
in the power spectrum shown in Fig. Hjo are sharper in the time direction than those derived 
in Fig. inb, thus indicating a better time resolution in the case of the Paul wavelet function. 
On the other hand, since the Morlet mother wavelet is narrower in the periodicity space than 
the Paul function (Fig. [7]), the same peaks appear less elongated and better defined in the 
period direction, thus indicating that a better frequency resolution is obtained if the Morlet 
function is applied (Figs. |3]d and[2}D). 

The wavelet transforms based on the Paul mother wavelet (characterized by few zero 
crossings. Fig. [7]) can in principle identify many types of signatures in the power spectrum, 
such as sharp variations in the input data like simple rounded impulses, and not just waves. 
Nevertheless, by taking into account that the same significant high-power features are found 
in both the Paul and Morlet analyses and that the Morlet function is more suitable to detect 
signatures of periodic and quasi-periodic fluctuations (since it contains more "waves" than 
the Paul function. Fig. [7]), it is possible to infer that the signals detected by CORl-A might 
be not simply variable. 

Thus, since the main aim of this study is to localize in time and in space any quasi- 
periodic feature associated with expected non-stationary processes, the Paul mother wavelet, 
which assures a better time localization than the Morlet function, has been used throughout 
the analysis reported hereinafter. 

In summary. Figs. [2)D,c,d show that significant non-stationary variability exists over 
a broad range of periodicities, extending from 0.5 up to 8 h; the strongest fiuctuations are 
found in the 4 — 8h band. The observed variability persists throughout the observational 
interval. The temporal evolution of the wavelet power is strictly related to the periodicities 
involved: longer periodicities are characterized by a more slowly evolving wavelet power. 

3.1. 2D reconstruction of the wavelet power 

The spatial and temporal evolution of the total brightness fluctuation power is obtained 
by computing the averages (Eq. lASp as functions of time, pixel by pixel, for the four selected 
periodicity bands. The results have been used to construct four different movies (available 
as online material in the electronic version of this paper), which show non-stationary, highly 
unstable fluctuating features in the entire FOV observed by STEREO- A CORl, appearing 
as arch-shaped (at low and mid latitudes) or ray-like (at higher latitudes) structures. The 
lifetimes of such features are comparable to the characteristic periodicities of each band. A 
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single frame of each movie at a given instant is derived by combining into a two-dimensional 
contour plot the values of the average wavelet power. Figure H] shows the 2D maps of the 
wavelet power in the four different periodicity bands, at the same instant as in Fig. [T| the 
extrapolated magnetic field lines (April 19, 2008 at 00:04 UT) are superimposed to the 2D 
maps. The figure also shows the SOHO/MDI disk magnetogram (as seen from the STEREO- 
A satellite) at that epoch. The power is plotted when exceeding the 99% confidence level 
(derived as specified in the appendix). The features identifiable in the different bands are 
significant non- stationary fluctuating structures, with power well above the background level. 
The background is mainly due to random fluctuations of the total brightness of the quiet 
solar corona. 

The features corresponding to power enhancements are often but not always localized 
along the magnetic field lines (Fig. H]). Moreover, they do not correspond to bright resolved 
structures with sharp edges in the total brightness images of the solar corona. The power of 
the bright features, shown in the Fig. HI is about ten times larger than that of the dimmer 
regions (Fig. |2]i). For instance, the power associated with the bright polar rays evident in 
the 4 — 8h band is about 2 — 4 x 10~^^ {SB/Bq)"^, whilst that of the dim regions between 
bright rays is 1 — 2 x 10~'^'^ {6B/Bq)^. Similar ratios are found for the wavelet power in 
the bright equatorial arches compared to that derived for the dark cavities. The periods of 
the embedded fluctuations are related to the thickness {d±) and the length (L) of the arch- 
shaped and ray-like features, determining the spatial scales of the fluctuating structures; 
that is, periodicities are longer for larger spatial scales. As an example, in the 4 — 8h band, 
the thickness of the fluctuating structures is d±_ ~ 0.01 -R© = 7 Mm, while the extrapolated 
lengths (down to the solar disk) are L ~ 2.9 Rq = 2.03 x 10^ Mm and L ~ 1.4 Rq = 980 Mm 
for the arches and rays, respectively. 



3.2. Spatial and temporal profiles of the fluctuation power 

In order to investigate in detail how the power varies and evolves along the bright 
non-stationary fluctuating features, spatial and temporal profiles of the power of the total 
brightness are derived as a function of a linear coordinate along arches and rays (Fig. H]). 
Since these structures are more evident for periods of 4 — 8 h, hereafter we focus our attention 
on this periodicity band. 

The left panel of Fig. |5] shows the power averaged over the 4 — 8h band, in the western 
hemisphere of the CORl-A FOV on April 19, 2008 at 09:35 UT. There is clear evidence for 
a system of two encapsulated closed arches, protruding from the occulted coronal region, i.e. 
from 1.4 .Rq, up to about 1.7 Rq and 2.0 i?©, respectively. The arches are separated by a 
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Fig. 4. — 2D maps of the total brightness fluctuation power averaged over the 0.5 — Ih, 
1 — 2h, 2 — Ah and 4 — 8h periodicity bands respectively, obtained in the FOV of CORl-A 
on April 18, 2008 at 21:50 UT; the extrapolated magnetic field lines and the SOHO/MDI disk 
magnetogram, as seen from STEREO-A on April 19, 2008 at 00:04 UT, are also shown; the 
plotted average variance of the total brightness is significantly above the 99% confidence for 
the assumed white noise; the dashed white line indicates the edge of the STEREO occulter. 
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dark cavity. The dots outline the outer arch. 

The spatial profiles of the fiuctuation amphtude along the outer arch are plotted as a 3D 
histogram in the top right panel of Fig. |3 The x-axis represents the curvilinear coordinate 
defined as the incremental distance along the path from the southernmost visible part of the 
arch. The y-axis represents the time (in hours) starting on April 18, 2008 at 16:35 UT. The 
spatial behavior of the total brightness fiuctuations along the arch is measured at each point 
by the percentage variation of their amplitude relative to the instantaneous total brightness 
(z-axis). Spatial profiles (along the x-axis) are obtained at a cadence of two hours in the 
period from April 18, 2008 at 16:35 UT to April 19, 2008 at 17:35 UT. The relative amplitude 
of the total brightness fiuctuations increases along the arch, reaching its maximum of about 
3% around the cusp of the closed structure, a factor ~ 2 higher than at the base. The trend 
of the spatial profiles is quasi-independent of the time. 

The temporal profiles of the fiuctuation power inferred at each point along the arch are 
plotted as a 3D histogram in the middle right panel of the same figure. The time interval 
over which the temporal profiles have been derived is much longer (more than 24 h) than the 
lifetime of the fiuctuating feature, thus assuring that its entire temporal evolution is indeed 
monitored. The x-axis represents the time (in hours) starting on April 18, 2008 at 16:35 UT. 
The y-axis represents the curvilinear coordinate along the arch. The z-axis is the normalized 
wavelet power for each temporal profile. The total brightness fiuctuations along the arch are 
typically characterized by a monotonic damping. This suggest that the fiuctuating structures 
found in the corona are non-stationary in nature. The 1 / e- width of the main fiuctuation at 
each point of the arch provides information on the lifetime of the fiuctuating structure, which 
is slightly less than ten hours, that is, of the same order of the temporal scales examined in 
the band of 4 — 8/i. The power reaches its maximum on April 19, 2008 at about 07:35 UT 
almost simultaneously along the arch, except at the northernmost base where the maximum 
is reached slightly earlier. This fact could be due to a tilt of the fiuctuating arch relative 
to the POS, which would produce a variable column depth along the LOS, thus modulating 
the temporal profiles of the wavelet power along the structure. 

The bottom panel of Fig. Oshows the temporal and spatial coordinates of the power peak 
in the oscillating arch. The error bars are derived from the 1 / e- width of the corresponding 
temporal profiles. The linear coefficient of the fit shown as a dotted line, B = 2.40 ± 
3.01 h Rq^, is consistent, within the uncertainties, with a fiat curve. This implies that the 
total brightness fiuctuation detected reaches its maximum power fairly simultaneously along 
the arch, hence the fiuctuations are coherent in phase. 

The same analysis has been also performed on a ray-like structure in the south-western 
hemisphere, shown in the left panel of Fig. |6l which shows the wavelet power averaged over 
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Fig. 5. — Left: Wavelet power averaged over the 4 — 8 /i periodicity band, inferred in the 
western hemisphere of the CORl-A FOV on April 19, 2008 at 09:35 UT; the dots outline 
an arch-shaped fluctuating structure; the continuous white line indicates the solar disk limb, 
whilst the dashed one the edge of the STEREO occulter. Top right: 3D histogram showing 
the spatial profiles of the variability amplitude relative to the instantaneous total brightness 
referring to different points along the arch, at different times (one every about two hours), 
from April 18, 2008 at 16:35 UT to April 19, 2008 at 17:35 UT. Middle right: 3D histogram 
showing the temporal profiles of the power of the fiuctuations inferred along the arch, in the 
same time interval as in the top right panel. Bottom right: Time instants when the wavelet 
power of the different temporal profiles reaches its maximum as a function of the curvilinear 
coordinate along the arch-shaped structure; the linear fit is shown as a dotted line. 
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the 4 - 8 /i band, on April 20, 2008 at 13:35 UT. 

The spatial and temporal profiles of the wavelet power of the total brightness along the 
ray are shown as 3D histograms on the top and middle right panels of the figure. They are 
very similar to those obtained for the arch-like structure and can be summarized as follows: 
i) the fluctuations strengthen with distance along the ray: within 1 Rq their amplitude, 
relative to the total brightness, increases by a factor of 2 up to 3%; further out the structure 
is no longer observed; ii) the fluctuations are non-stationary and damped with a lifetime of 
about ten hours; the lifetime of the structure is of the same order of the periodicities of the 
embedded fluctuations, that is slightly less than ten hours; iii) the fluctuations along the 
ray- like structure are coherent in phase. 



4. Discussion 

Strong non-stationary fluctuating structures have been detected in the solar corona, as 
closed arches in the equatorial regions and ray-like features in the polar holes. They are often 
found in correspondence with the strongest lines of the coronal magnetic field as shown in 
FigHl The observed variability can be ascribed to variations in the electron density of the 
emitting coronal plasma, since the white-light emission in corona is due to the Thomson 
scattering of the solar disk photons by the free electrons. The periodicities observed in the 
fluctuating corona cover a wide range of scales, from 0.5 to 8h; the strongest signals are 
found in the 4 — 8h band, where the fluctuating structures have a characteristic transverse 
and longitudinal spatial scale of ~ 0.01 Rq = 7 Mm and ~ 1.4-2.9 Rq = 980-2030 Mm, for 
the ray and arch features, respectively. These resu lts confirm the findings of th e first studies 



on dens ity fluctuations in the outer solar corona by lBemporad et al.l (120081 ) and lTelloni et al. 



(l2009al Jbl). In these papers, quasi-periodic oscillations with periodicities larger than a few 
hours were found in the Fourier power spectra; however, the analyses were much more limited 
in space and time, being performed with the data obtained with a spectrometer. 

Notwithstanding we have ascertained that the total brightness variability is non-stationary 
in nature, it reveals a long-term persistence throughout the observational time interval of 9 
days. This fact suggests that a non-transient, persistent mechanism of quasi-periodic nature, 
which not necessarily acts at coronal level, continuously drives density fluctuations in the 
solar corona. The fluctuations continuously originate and successively decay in a time that 
is of the same order of their characteristic temporal scales; this causes a modulation of the 
variance of the total brightness as seen in Figs. |2^ and [2jl and in the maps of the corona 
of Fig. HJ Such a behavior is clearly evident in the movies of the spatiotemporal evolution 
of the wavelet power (electronic version of this paper). It can be suggested that the density 
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Fig. 6. — As in Fig. [5] for a ray-like fluctuating structure (outlined by dots) in the south- 
western quadrant, on April 20, 2008 at 13:35 UT. The spatial and temporal profiles refer to 
the time interval from 05:15 to 17:20 UT on April 20, 2008, sampled every about one hour. 
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fluctuations, which are initially distinct and characterized by a single dominant periodicity, 
either rapidly decay into a turbulent mixed flow via nonlinear interactions with other plasma 
modes, or are damped by thermal conduction. 

Furthermore, the results obtained in this study point to a tight correlation of the tem- 
poral evolution of the power of the plasma variability detected in this work and the spatial 
and temporal scales of the fluctuating coronal features. Longer periodicities are indeed typ- 
ical of fluctuating structures of larger size and evolving more slowly. This reflects an overall 
characteristics of the Sun: larger-scale solar structures have characteristic timescales longer 
than those corresponding to smaller-scale structures. For instance, the strong correlati on 



existing between the size and the time duration of the polar plumes (iLlebaria et al.ll2002[ ) is 
well known. This behavior is as well evident in the movies (electronic version of this paper) 
that show the spatial and temporal evolution of the magnetic features, both arches and rays, 
in the corona from 1.4 to 4.0 Rq. 

Finally, the analysis of the temporal evolution of the wavelet power along the non- 
stationary arch-shaped and ray-like structures (§ 3.2) reveals that the plasma fluctuations 
are highly coherent, since each point along the coronal magnetic structure fluctuates with 
the same phase. The phase coherence of the plasma fluctuations might indicate that they 
are likely to be generated simultaneously by the same driver mechanism. 



4.1. Assessment of solar rotation effects 

Even if the frequency contribution due to the solar rotation has been filtered out, in 
order to assure that the same volume of the corona is analyzed (as described in § 3), it cannot 
be a priori excluded that, in observations performed at the limb, tangential discontinuities 
separating adjacent radial structures with different electron density (for instance, magnetic 
flux tubes CO- rotating with the Sun), could mimic wave phenomena. During limb observa- 
tions of the quiet corona, the maximum density along the line of sight is at the plane of the 
sky, due to its proximity to the disk. Thus a density variation caused by quasi-regularly 
spaced tangential discontinuities, due to the rotation of the Sun, induces a quasi-periodic 
signal in the coronal emission detectable at best when crossing the plane of the sky. We 
have tested this occurrence by evaluating the temporal scale of rotating long-lived struc- 
tures, Tr = d±j^;^, where d± is the geometric thickness of the oscillating arch, perpendicular 
to the magnetic field lines, Tr ~ 27.5 d is the rotation period observed in the outer corona 



( IGiordano fc Mancusdl2008l ) and r ~ 1.5 Rq is the average heliocentric distance of the fluc- 
tuating coronal arch. By assuming d± ~ 0.01 Rq, that is the measured value in the 4 — 8h 
periodicity band, we obtain < 1 /i, which is signiflcantly smaller than the periodicities 
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characterizing the fluctuating arch. We have also addressed this issue by assessing the effect 
of the rotation on the temporal evolution of a c oronal magnetic loop, formed by numerous 
spatially unresolved filamentary sub-structures (iKlimchukl 120091 ) . The loop fine structure 
has been modeled by considering a curved cylinder of longitudinal and transverse dimen- 
sions of 2.9 i?0 and 0.01 i?0, respectively (dimensions consistent with the length, L, and the 
thickness, dj_, of the observed arch-shaped structure), consisting of many sub-elements char- 
acterized by different densities. On the basis of this model we compute the coronal electron 
density as a function of time at the plane of the sky in a time interval of about one day. In 
the hypothesis that the randomly distributed sub-structures produce a correlated signal, the 
Fourier analysis of the synthetic coronal emission shows that the dominant mode of variabil- 
ity is consistent with the crossing timescale of the rotating arcade-system, i.e. ^ Ih, whilst 
significantly less power is found at longer periods. 

The effect of the solar rot ation should be neg ligible in the case of polar structures 
which develop at high latitudes (jWilhelm et al.ll201ll ). For instance, the s hape of the plume s 
on spatial sca l es of ~ 30 arcsec is nearly invariant for hours to days (jWaldmeieii Il955l ). 



Llebaria et al.l (120021 1201 ll ) describe polar plumes as enduring structures showing brightness 
variations which might be ascribed to the solar rotation at timescales longer than about 10 h. 
The intrinsic variations of the plume brightness on scales of the order of 6 h, due probably 
to the intermittent nature of the generating process, are superimposed on this long-term 
variability. Therefore, in the 4 — 8h band, where, in a plume-like feature, we find significant 
variability, the effect of solar rotation should be negligible. 

Before concluding that the arch-shaped power enhancements can be interpreted as due 
to intrinsic coronal density fluctuations or even more MHD waves, we have also to rule out 
another effect that could in principle mimic power in the wavelet maps, i.e. any possible 
spatial motion of an arch-shaped structure with clear edges. Quasi-periodic spatial displace- 
ments of coronal loops, detected by TRA CE in the inner corona, have been reported in 
many papers (e.g. lAschwanden et al.lll999[ ). They are commonly triggered by nearby flare 



events and are typically kink instabilities characterized by timescales of the order of a few 
minutes. Note, in this respect, that no flares or eruptive phenomena were detected during 
the examined observing run and, moreover, that the fluctuations found in this study are 
characterized by much longer timescales. Finally, power due to non-intrinsic variability such 
as that generated by a simple pulse of material injected in the arch with a duration appro- 
priate to trigger the timescales of the detected fluctuations can be also ruled out on the basis 
of the high coherence of the signal found all over the characteristic fluctuation lifetimes. 



The analysis shows that the non-stationary density fluctuations outlining closed and 
open coronal structures, such as arch-systems and ray-like features, can be ascribed to in- 
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trinsic brightness variations with possible MHD wave character. The compressive compo- 
nents of MHD waves can indeed perturb the ambient solar corona, modulate periodically or 
quasi-periodically the electron density and, in the end, generate plasma fluctuations, which 
can be observed by a white-light imaging coronagraph such as STEREO- A CORl. 



5. Interpretation 
5.1. Arch-shaped structure fluctuations 

Theoretically various MHD waves are pr edicted in coronal loops. The dispersion relation 



for magneto-acoustic waves (e.g. lPriestlll987l ) is satisfied by two kinds of waves: a slow-mode 
branch, with acoustic phase speeds, and a fast-mode branch, with Alfven speeds. The class 
of the fast magneto-acoustic waves includes a symmetric sausage mode and an asymmetric 
kink mode. In the sausage mode, opposite sides of an arcade oscillate in anti-phase and 
the central axis remains undisturbed, while in the kink mode, the opposite sides and the 
central axis oscillate in phase. However, both modes drive plasma fluctuations characterized 
by the same phase along the arcade, as those observed in the coronal arch and detected in 
the present analysis (Fig. Hence, these MHD modes, as well as slow-mode waves, can be 
considered physical candidates for the interpretation of the observed variability. 



5.1.1. Fast-mode waves 



Kink modes produce transverse amplitude oscillations. Oscillatory transverse displace- 
ments of coro nal loops have been observed i r i the inner corona in EUV wavelengths with 



TRA CE (e.g. lAschwanden et al.l Il999l . l2002l : iNakariakov et al.l Il999l : ISchriiver et al.l Il999 



2OO2I ) . In the hypothesis that a standing kink-mode is responsible for the observed electron 
density fluctuations, the characteri stic period P k ink of such a mode at the fundamental har- 
monic can be derived, according to lAschwandenI (120041 ). as Pkink ~ ~ where and va 
are the kink and the Alfven speed, respectively. In low-/3 plasma conditions, where the ther- 
mal pressure is much smaller than the magnetic pressure, ^ va, with va ~ 1000 kms~^, 
for typical electron densities and magnetic fields in the coronal region of interest. For the 
observed arch with a measured length L = 2.03 x 10'^ Mm, the characteristic period is 



kink 



Ih. 



Hence, if the fluctuations detected in the coronal arch were due to standing kink waves, 
we would expect kink-mode periods, and in turn electron density fluctuations, of less than 
one hour. It turns out that, according to the analysis, the observed variability cannot be 
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ascribed to a standing wave in the kink mode with fixed nodes at the footpoints. 



Sausage MHD modes can be excluded as well because they are expected to be completely 
suppressed at small wavenumbers. In a coronal arch with length and cross section of the 
order of the observed ones (L ~ 2.03 x 10^ Mm and d± = 0.01 i?©,) the dimensionless 
wavenumber is of the order of 0.005, much smaller t han the cutoff wav enumber for typical 
coronal conditions, which falls into the range 0.8 - 2.4 (jAschwanderul2004j ) . As a consequence, 
the sausage type waves cannot explain the density fiuctuations observed in the arch-shaped 
structure. 



5.1.2. Slow-mode waves 

Slow-mode oscillations in coronal loops have been imaged for the first tir ne in the inner 
corona with the SUMER (Solar Ultraviolet Measurement of Emitted Radiation. IWilhelm et al, 
19951 ) spectrom eter on SOHO, by measu ring the Doppler shift and intensity variations of 



EUV lines (e.g. IWang et al.ll2002l |2003aJJb|). Slow-mode acoustic waves are essentially com- 
pressional, longitudinal sound waves. In coronal conditions, for typical temperatures of 
T ~ 1.5 Mi^', the sound speed Cs ~ 180 km s"^ is generally much smaller than the Alfven 
speed {va ~ 1000 km s~^), that is Cs/va ~ 0.2. Thus, the phase speed of slow magneto- 
acoustic waves in coronal loops, namely the tube speed ct, 



Ct 



CsVa 



for Cs <^ va 



(1) 



is close to the sound speed Cg (e.g. |Priestlll987( ). 



The period of a MHD slow-mode standing wave, having the two endpoints of a magnetic 
fiux tube as fixed nodes, is the sound travel time of the fundamental harmonic, back and 
forth in the coronal loop, Pgiow = ~ 7^- Hence, the period of a slow magneto-acoustic 
standing wave, Psiow ~ 6 /i, supported by the observed fiuctuating arch structure falls in 
the 4 — 8 /i band where most of the power of the detected fiuctuations is observed; thus, 
the fiuctuation of the magnetic arch extending in the outer corona out to approximately 
2 Rq is consistent with a slow magneto-acoustic standing wave. Therefore, the MHD slow- 
mode standing waves are ideal candidates for the interpretation of the plasma fiuctuations 
embedded in the coronal arch detected in the equatorial region of the quiet corona. 

As far as the origin of the slow-mode wave is concerned, one possibility is that the vari- 
ability is triggered by a pressure disturbance locahzed near one of the footpoints, which prop- 
agates as a slow-mode magneto-acoustic wave in longitudinal direction along the arch and 
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is reflected at the opposite footpoint, thus generating a standing wave (a compression back 
and forth along the magnetic flux tube). A stochastic excitation by convective supergranular 
motions generates slow magneto-acoustic waves with a period which is consistent with those 
observed in the coronal fluctuating arch. The lifetimes of supergranular cells are ranging 



from 6 to 20 /i (lAigrain et al.l l2004j ) : this interval is partially overlapping with the A — 8h 
band, where closed magnetic arches clearly exhibit the largest wavelet power. A further sup- 
port to this interpretation is given by the value estimated for the wavelengths of slow-mode 
standing waves triggered by supergranulation, A,, = ctTs ~ CgTs ~ 5.6 Rq = 3.89 x 10^ Mm. 
This wavelength is approximately twice the loop length L of the fluctuating coronal arch, 
hence, it is consistent with a standing wave with fixed nodes at the footpoints of the mag- 
netic flux tube. The present discussion suggests that supergranular convective motions very 
likely play a crucial role in the excitation of the slow MHD standing waves trapped in the 
large closed fluctuating magnetic structures extending in the outer corona. The length of 
the arch studied in this analysis is compatible with a spatial scale that can trap the standing 
wave, thus acting as a filter of standing waves with the appropriate wavenumber. 

The suggestion that the driver of the standing waves along coronal arches can be identi 



tied wit h the supergranular convective motions is consistent with the results by lTelloni et al. 



(l2009al ). where the characteristics of the coronal density fluctuations observed at 1.7 i?© in 
low-latitude regions, such as the temporal and spatial scales, the degree of persistence and 
phase coherence, were all indicating a likely correlation between the variability of the so- 
lar corona and the p hotospheric dynamics induced by supergranular convection. However, 
Telloni et al.l (l2009al ) interpret the correlation between the dynamics of supergranulation and 
the coronal variability by assuming that structures of the supergranulation at the lower levels 
of the solar atmosphere is maintained in the outer corona; thus, coronal structures, crossing 
the POS during rotation, produce a periodic signal, whose characteristics depend on the 
temporal and spatial scales of the supergranulation. In this study, the link to convective mo- 
tions is different, since, according to the interpretation of the present results, the stochastic 
motion of the supergranular cells acts as a driver for the excitation of slow magneto-acoustic 
waves which remain trapped in the arch structure as standing waves when their wavelengths 
are compatible with the loop length. 



5.1.3. Damping of MHD modes 

The MHD slow-mode standing waves, inferred from the density fluctuations of the arch 
observed in the quiet corona, show a strong damping effect, which has a decay time of 
the order of the fluctuation period, as shown in the middle right panel of Fig. [51 This 
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rapid damping can be ascribed to processes such as thermal conduction, radiative losses 
and phase mixing. Thermal conduction and radiative losses are the major non-adiabatic 
cooling effects, which can contribute to the damping of the wave motion. For typical coronal 
electron densities (rie ~ 1.2 x 10^ cm~^) and temperatures (T = 1.5 MK), at the average 
heliocentric distance (r ~ 1.5 Rq) of the fluctuating arch of length L ~ 2.03 x 10^ Mm, 
the cooling time due to thermal con duction from the corona to the chromosphere is of the 



order of Tcond ~ 5 — 6 /i ( jPriestl 119871 ). Therefore, thermal conduction can be considered an 
important factor in damping standing slow magneto-acoustic waves trapped in the closed 
fluctuating coronal arch of Fig. [5l since the conductive cooling time is of the same order of the 
wave damping time. A — 8h. The interpretation of thermal conduction as the main physical 
mechanism responsible for damping standing slow magne to-acoustic waves in ma gnetic loops 



is also supported by the MHD simulations performed by lOfman &: Wangl (j2002[ ) in order to 



interpret SUMER observations. Whilst, since the radiative cooling time, defined by th e ratio 



of the thermal energy and the radiative cooling rate, is generally lower than one hour (jPriest 



19871 ). radiative losses are negligible in damping the MHD slow- modes, which have larger 



decay times. 

Another mechanism that could substantially damp MHD slow-mode waves is phase 
mixing. Magneto-acoustic waves in a homogeneous medium, in density and temperature, 
propagate undamped and undisturbed because they are in perfect resonance with each other. 
However, coronal loops have large density and temperature gradients, especially near the 
footpoints where the standing waves are reflected. For instance, for MHD slow-mode waves, a 
density decrease along the coronal loop would increase the sound speed and decrease the wave 
periods along the magnetic flux tubes. Therefore, standing slow magneto-acoustic waves in 
a coronal loop might suffer intense phase mixing, during which modes of same frequency 
travel with different speed and then have different wavenumbers. As they propagate they 
become quickly mutu ally out of phase , leading to a destructive interference, and thus to 



wave dissipation (e.g. iHood et al.l 119971 ). In other words, the slow MHD modes, which are 
initially distinct and of a single dominant periodicity, might rapidly decay into a turbulent 
mixed flow via nonlinear interactions with other reflected modes. 

In conclusion both thermal conduction and phase mixing are viable mechanisms to 
explain the damping of the observed standing waves which manifest themselves in closed 
magnetic structures. 
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5.2. Ray-like structure fluctuations 



In analogy to the proposed interpretation of the electron density fluctuations in coronal 
arches, we have evaluated the possibility that outward propagating slow MHD waves might 
be responsible as well for the non-stationary plasma fluctuations detected in a polar region 
along the magnetic plume-like feature shown in the left panel of Fig. El which might have 
an open magnetic configuration. 

The density modulation due to propagating MHD waves would result in time-shifts of 
the temporal profiles of the wavelet power inferred along the polar ray. This is because nodes 
propagate with the wave. Hence, the density fluctuations are expected to be out of phase 
at different heliocentric altitudes along the open structure which might act as a channel for 
outward propagating waves. However, different part of a ray-like structure are observed to 
fluctuate in phase, reaching their maximum power quite simultaneously (middle and bottom 
right panels of Fig. E]). No evidence for magneto-acoustic wave propagation in the polar 
ray-like feat ure has been be fou nd by studying the temporal cross-correlations (method 
described by iDuvall et al.l ( 1l993l )) of the total brightness fluctuations observed at different 
locations along the structure. This clearly indicates that the plasma density fluctuations 
detected along the ray of Fig. [6] have to be ascribed either to quasi-periodic motion forth 
and back across the POS of the entire structure or to intrinsic brightness variation of the 
polar ray itself, likely due to the intermittent nature of the process underlying its generation. 
As a matter of fact, although the magnetic configuration of solar polar regions is relative 
stable, the mechani sm driving t he formation of the actual ray structure might operate only 



intermittently (e.g. IWanglll998[ ) 



Based on the characteristic timescales of the fluctuations detected in the present analysis, 
the main mechanism responsible for the temporal variations observed in the brightness of 
the polar features might be caused by magnetic reconnection which should be acting with 
comparable timescales. It has been suggested that magnetic reconnection of continuously 
emerging flux, guided by supergranular convection, could affect significantly the li fecycle 
of polar plumes, driving multiple hourly occurrences of the plume brightening (e.g. IWang 
19981 : iRaouafi et al.ll2008l ). Magnetic reconnection in polar regions has also bee n suggested 



as source of waves and energy to explain the fast wind (e.g. lAxford et al.lll999[ ). The polar 
rays observed in coronal maps of the wavelet power not necessarily are observed in coronal 
brightness maps, thus they might not always coincide with polar plumes and they are more 
extended than plumes. The same process of magnetic reconnection can in principle act at the 
base of plumes as well as rays, when they do not coincide. Hence, the lifetime of the density 
fluctuations derived in the present study (of the order of ten hours, middle right panel of 
Fig. E]) can provide information both on the duration of the intermittent phenomena and 
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on the lifecycle of the plume itself, provided the same process is acting at the base of both 
kinds of structures. Moreover, the characteristic lifetime of the detected fluctuations could 
also be related to the timescal e corresponding to the pronounced teri dency of recurrence of 
plumes at the same locations (ILamy et al.l 119971 : iDeForest et al.ll200ll ). 



6. Conclusions 

This paper reports on the detection of transient periodic or quasi-periodic phenomena 
occurring in closed and open magnetic structures of the quiet solar corona out to 2 Rq. The 
technique used to reveal and characterize the fluctuations is based on the wavelet trans- 
forms, which represents a powerful tool to investigate the fine spatial pattern and temporal 
variability of fluctuating coronal features. Here, this technique has been used to map the 
variance of the total brightness of the outer solar corona from 1.4 to 4.0 Rq for the first time 
and in a systematic way. 

The results show that the solar corona is permeated by quasi-periodic recurrent bright- 
ness fluctuations, characterized by a wide range of spatial and temporal scales. They are 
however not uniformly distributed but clearly outline the coronal magnetic field, thus ap- 
pearing as embedded in fluctuating magnetic structures. Moreover, they are significantly 
enhanced in regions of stronger field lines. In particular, it has been possible to identify 
an extended coronal arch in the mid-latitude equatorial region, likely part of a complex 
equatorial streamer present at the West limb, and a ray structure at high latitudes, both 
characterized by fluctuations of large amplitude. The arch-like structure is not observed in 
total brightness maps, thus implying that it is less dense than adjacent coronal structures. 
In the closed low and mid-latitude magnetic arch, the observed plasma fluctuations are con- 
sistent with the presence of slow magneto-acoustic waves trapped as standing modes and 
damped on about ten hours. These waves are very likely excited in a stochastic way by 
the convective motions at supergranular scales. The most likely mechanisms which could 
explain the wave damping are thermal conduction and phase mixing. In the polar region, 
the density plasma fluctuations of the ray structure, not necessarily coincident with a plume 
and even more extended than a plume, might be interpreted as signatures of intermittent, 
quasi-periodic magnetic reconnection due to continuously emerging flux tubes. This mech- 
anism, considered the most probable one for the generation of plumes, as suggested in the 
literature, can in fact induce waves and ejections of mass; however, in order to explain the 
observed fluctuating feature, it should occur in a quasi-periodic recurrent way. 

In conclusion, the present study of coronal density fluctuations, which extends previous 
analyses on periodic or quasi-periodic phenomena to cover the entire solar corona from the 
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limb out to 4i?0, provides new further information on the MHD coronal wave properties. 
The outer solar corona is a region where wave observations have been up to now highly 
fragmentary in space and time. The measurements of the power of MHD waves in the outer 
solar corona performed in this paper provide an accurate spatiotemporal reconstruction of 
quasi-periodic plasma phenomena. These may lead to a significant progress in determining 
how the energy carried by MHD waves is transferred from the photosphere to the corona, 
and there dissipated to accelerate the solar wind. 
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The authors are thankful to the anonymous referee for helpful suggestions leading to a 
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A. Wavelet transforms and wavelet power spectrum in the analysis of the 

solar data 

A.l. General formalism 

Given a wavelet function, ^0(^7)) depending on a non-dimensional 'time' parameter 
77, that has zero mean and is localized in both time and frequency space, the continuous 
wavelet transform of a A'"-point discrete two-dimensional time series /(r, 0, t), sampled at 
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equal time spacing 5t, where r and cj) are for instance the spatial polar coordinates and 
t = 0, (N — l)St is the time, is t he convolution of / with a scaled and translated version 



of "^oiv) (ITorrence &: Compdll998l ): 



{N-l)5t I 

mr,0,t,.) = -^ /(^''^''^')(7)^o(^)> (Al) 

^ t'=0 

where the asterisk indicates the complex conjugate and ^'o('7) is normalized to unit 
energy. By varying the wavelet scale parameter s and translating the wavelet function 
\E'o('7) in the time domain, it is possible to determine the amplitude of any periodic feature 
contained in the data and how this amplitude varies with scale and time. It is worth noting 
that the wavelet scale s does not necessarily correspond to the Fourier period, r, because 
of the functional form of the wavelet function. The relationship between the equivalent 
Fourier period and the wavelet scale is derived analytically for a particular wavelet function 
by substituting a cosine wave of a known frequency u = into Eq. I All and computing the 
scale s at which the wavelet power spectrum reaches its maximum. 

There are a number of different wavelet functions ^E'o(^) with different properties that 
can be best suited to be applied to different studies. The best choice depends on first instance 
on its shape and width and whether it is a complex or real valued function. Complex wavelets 
are ideal candidates for capturing oscillatory behavior with respect to real wavelets, which 
are more useful to isolate peaks and discontinuities. The width of the wavelet function, 
namely the e-folding time of the wavelet power, is related to the resolution of the wavelet 
itself: a narrow function in time has good time resolution but poor frequency resolution and, 
viceversa, a broad function has poor time resolution and good frequency resolution. Two 
common mother wavelet functions, both complex, are the Paul and Morlet wavelets. The 
Paul wavelet is defined as 

MV) = ^T^il - ^Vp-^'\ (A2) 

where m (= 4 in the present analysis) is the order of the wavelet. The width is s/\/2 
and the wavelet scale s is related to the Fourier period r by the expression r = 2m+i ^ ~ lAs. 

The Morlet wavelet consists instead of a plane wave modulated by a Gaussian: 

^o(ry) = 7r-i/^e'"o''e-'''/2^ (A3) 
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where uq (= 6 in the present analysis) is the nondimensional frequency. The width 
is a/2s and the wavelet scale s is related to the Fourier period r by the expression r = 
^^^s = 1.03s. 

The Paul mother wavelet is more sharply defined in time in comparison to the Morlet 
function. Moreover it is characterized by few zero crossings and for this reason it can in 
principle identify many types of signatures in the power spectrum, as sharp variations in the 
input data and not just oscillatory modes. The Morlet function is narrower in periodicity 
space than the Paul wavelet, leading to a better frequency resolution. Moreover it contains 
more oscillations and then it is optimized for detecting oscillatory signals. Pictures of the 
real (solid lines) and imaginary (dashed lines) parts of the Paul and Morlet wavelets in both 
the time domain (top panels) and the frequency domain (bottom panels) are shown in Fig. 

m 

The wavelet power spectrum is defined as the square of the modulus of the wavelet 
transform \W{r,(f),t, s)\'^. A vertical slice through the wavelet power spectrum is a measure 
of the spectrum at a given instant, so that the average of the wavelet power spectrum over 
the whole observational period, < t < {N — l)6t, gives the global wavelet spectrum, defined 

as 

|W(r,0,s)p = - J2 \W{r,^,t,s)\'. (A4) 
t=o 

It is worth noting that the global wavelet spectrum appro aches the resul ts that are 



obtained with the Fourier spectrum of the considered time series (lPercivallll995l ) 



The average wavelet power, that is the sum of the power over all the considered peri- 
odicities between si and S2, is given by: 



where 6j (=0.125 in the present analysis) is the spacing between the discrete scales 
commonly defined as fractional powers of two, Sj = sq2^^^ (sq is the Nyquist- Shannon scale), 
and Cs is the reconstruction factor for the assumed wavelet function (we adopted 1.132 and 
0.776 for the Paul and Morlet wavelets, respectively). 

The scale-averaged wavelet power defined by Eq. IA5I gives a measure of the average 
variance of the time series as a function of time, over the range of periodicities considered in 
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Paul wavelet basis Morlct wavelet basis 




sgj/(27t) scj/(27t) 

Fig. 7. — Top left: Real part (solid) and imaginary part (dashed) of the Paul wavelet (m = 4) 
in the time domain. Bottom left: Paul wavelet in the frequency (a;) domain. For plotting 
purposes, the scale s was chosen to be s = 2505t. Top right: Real part (sohd) and imaginary 
part (dashed) of the Morlct wavelet (cuq = 6) in the time domain. Bottom right: Morlet 
wavelet in the frequency (a;) domain. For plotting purposes, the scale s was chosen to be 
s = 2bQ5t in both panels. 
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the analysis. 

A. 2. Cone of influence 

Because of the finite length of the time series, errors will arise at the beginning and end 
of the wavelet power spectrum, as the Fourier transform assumes the data are cyclic. The 
Cone Of Influence (CGI) is deflned as the region of the wavelet spectrum where edge effects 
become important and is commonly defined as the e-folding time for the autocorrelation of 

the wavelet power at each scale. This c-folding time is chosen so that the wavelet power for 
discontinuity at the edge drops by a factor and assures that the edge effects are negligible 
beyond this point. 

A. 3. Significance levels 

Determining significance levels for a wavelet spectrum needs to choose an appropriate 
background spectrum in order to establish a null hypothesis against which the inferred spec- 
trum must be tested. If a peak in the wavelet power spectrum is significantly above the 
background spectrum, it can be assumed to be a true feature with a certain percent con- 
fidence. In many astrophysical applications it is a good choice to adopt as a background 
spectrum a fiat spectrum resulting from wavelet analysis of pure white-noise, normally dis- 
tributed time series. The wavelet coefficients of a normally distributed random variable are 
normally distributed too; the square of a normally distributed variable is chi-square dis- 
tributed with one degree of freedom, Xi- Hence, in the hypothesis of white- noise normally 
distributed time series, the wavelet power spectrum would be x| distributed with two de- 
grees of freedom. The 99% confidence contour lines are then determined by multiplying the 
background spectrum by the 99th percentile value of xh ^^^^ xli^^^^)- 



